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The method of steepest descents is extended to the case when ;i saddle point and a 
pole of arbitrary order are involved. An application to a problem in diffraction theory 
is demonstrated. 



1. Introduction 

The solution of numerous problems in electro- 
magnetic and acoustic wave propagation theory is 
frequently given by an integral representation of 
the form 



<p(z)= i w(v)e- 



]) dv. 



(1) 



Here <p(z) denotes the desired property (for instance 
the electromagnetic or acoustic potential). The in- 
tegration in (1) is taken along a curve C in the 
complex v plane along which w(v) and u(v) are 
defined. These function contain generally certain 
parameters which depend on the physical configura- 
tion under investigation. The parameter z is usually 
related to the distance between some point of refer- 
ence and the point of observation. A reduction of 
the integral (1) to an asymptotic series lias to be 
carried out when approximate expressions for large 
z are desired. The frequently used method of steep- 
est descents [3, p. 86] 1 becomes invalid when one of 
the above mentioned parameters assumes a value 
for which a pole of w(v) coincides with a saddle 
point of u(v). Methods to cope with such cases 
have been worked out and discussed [1], [8], [9]. 
The former two give expansions which proceed in 
terms of hypergeometric functions besides a "leading 
term", while the latter consists of an error integral 
as leading term and an asymptotic series with in- 
verse powers of z. It is shown here that the method 
given in [9], which involves a pole of the first order 
and a branch point of order 2 can be extended to 
the case involving a pole and a branch point of 
arbitrary order. The expansion for this case con- 
sists of a finite number (equal to the order of the 
pole) of terms involving Whittaker's functions and 
an asymptotic series with inverse powers of z. 

2. A Modification of Watson's Lemma 

It is assumed that the reduction of the contour 
integral (1) (for instance by the method of steepest 
descents) to one or several integrals of the Laplace 
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transform type 



m 



= ( F{t)e-'- 

Jo 



dt 



(2) 



has already been carried out and only this kind of 
integrals shall be considered. A theorem commonly 
referred to as Watson's lemma [2, p. 218] admits the 
asymptotic representation of f{z) for large z when 
the behavior of F(t) for small t is known. In its 
general form given by Doetsch [3, p. 48] and [4] one has 
the 

Theorem: Let F(t) be analytic in a sector a< 
arg t<CP of the complex t plane with the possible ex- 
ception of the origin, and let the behavior of F(t) in 
this sector be 

(a) \F(t)\<Ae c[t[ , for \t\>R, A, c>0, 

00 

(b) F(0«S^, -l<Rca <Rea 1 <- • ■ 



tohen t->0 in a< arg /<£. 



The asymptotic expansion for f(z) = F(t)e~ zt dt 

for large z is then obtained by inserting the series (b) 
for F(i) into (2) and integrating term by term. 

f(z) = i F(t)e- tz dt^f:a n Ta+ccn)z-- n - 1 

JO n = 



(3) 



when 2-^coin---— /3<arg 2<i 



The function F(t) in (2) is now assumed to be of 
the form 

F(t) = t«g(t), Bea>-1, (4) 

where the function g(t) is analytic in the vicinity of 
2 = 0. In this case 



with 



/w =i" 



t"(j(t)e-' z dt 



(5) 



g{t)=1LaJ.", a n =±-J£- n g(t)] . ( 6 ) 
97=0 ni \_jii j t=0 



15 



t 2 



to 



This series is convergent for |£|<3o|, where t is that 
singular point of g(t) which is closest to the origin. 
The other singular points are denoted by t u t 2 , . . . 
such that their moduli form a nondecreasing se- 
quence. We assume that none of the singular points 
is on the positive real £-axis. The asymptotic ex- 
pansion of (5) then becomes by (3) and (6) 



m 



-;; 



t a g{t)e- tz dt^ 



*ZXr(a+rH-i)2-«- 

71 = 



(7) 



for z— > oo . The range of validity for arg z depends on 
the behavior of g(t) according to (a) and (b) of the 
theorem quoted before. The case in which the 
singular point t is a pole of order m, say, is considered 
now. The Laurent expansion of g(t) in the vicinity 
of t is then 

g(t)=b- m (t-to)-»+ . . . +b- x (t-t )- l +g*(t) (8) 
where 



g*(t)=g(t)-b- m (t-to)- n - • • • -hS-toY 



(9) 



is analytic for |£|<C|£i| and consequently can be 
represented by the series 



«=0 



(10) 



It has to be remembered that in physical applications 
the function g(t) in (4) contains generally certain 
parameters which depend on the configuration of the 
problem, and this may become such that the pole 
t of g(t) coincides with the origin t=0 (see section 3). 
In this case all the coefficients a n tend to infinity by 
(6) and (8) and the asymptotic series in (7) becomes 
useless although the integral in (7) may exist. It is 
therefore necessary to obtain an expansion which for 
t ^0 represents the integral (7) asymptotically for 
large z and which for t =0 (if the integral exists) 
gives its correct value. Such an expansion can be 
obtained as follows: One inserts the representation 
(8) for g{t) into (7) and gets 



r* oo 

/(«)= J t a g{t)e- tl dt 

m f* °° n oo 

= S&-* (t-to)~H a e-' z dt+\ t a g*(t)e-"dt. 
*-l Jo Jo 



But 



(11) 



/; 



t«e- lz (t-t )~ k dt=T(l+a) (-to) l a -i k z ik ^ _1 



X e~ ^ W_ ik _ u . p+ »(- U z) , (12) 



by virtue of the integral representation for Whit tak- 
er's function [5, p. 274]. The integral 



r* oo 

Jo tC * g * 



(t)e- tz dt 



can be represented asymptotically for large z by 
means of the power series expansion (10) (which is 
valid for |*|<|*i|) and Watson's lemma. One obtains 
then finally 

m=J"t"g(t)e- tz dt 

^S^r(a+n+i)3-"- n - 1 +r(i+a)(-^ )^-^- 1 



72=0 



Xe~ z 



(13) 

The range of validity for arg z depends on the 
behavior of g*(t) according to the requirements (a) 
and (b) of the theorem quoted before. For certain 
combinations of a and m (13) can reduce to a simpler 
expression, as for instance 

W kti (z)=2l- k ziD 2k ^[(2z)i], 

(parabolic cylinder function) (14) 

W 7 _i,ife)=7rM^ 2 Erfc (s*) = zi«* f TV**" 1 * (15) 
W_i_ }0 (z) = -z-hi z Ei(-z) = z-h* z f t- l e- l dt (16) 

W ^(z) = (z/t)*kJ~z\ (modified Bessel function). 

(17) 

The case a= — \ ) m—\ (branch point of order 2 and 
pole t of the order 1 ) which leads to the special case 
(15) has been treated before [9]. For this case one 
obtains from (13) and (15) 



16 



f(z)=j Q "rig(t)dt ~ij c n T (n+g) *—* 

+ &_,*(- < n )-V- z 'oErfe[(-^ ) J]. (18) 

The coefficients c» according to (9) and (10) arc 
defined by 

09) 
(20) 



m-b-i(t-t )- i =T,c n t n , 

71 = 

or 

But these coefficients remain finite when t - 



>0. 



3. Application to the Wedge Diffraction 
Problem 

The investigations outlined before will now he 
applied to a previously discussed problem. It was 
shown [6, 7] that the solution of the problem of the 
diffraction of plane waves by a wedge or corner is 
governed by the integral 



ith 



I («,+)= f"rig(t)e- ypt dt 



•{cosh [Tr/a\0g(l + t+(t 2 + 2t)^)]-V0s(7nP/a)] 



(21) 



(22) 



Here g(t) is analytic near the origin and the singu- 
larity nearest to the origin is a pole of the order one. 
The pole t and its residue b ., are found to be 



* =-2sin 2 (#) 



(23) 

(24) 



The parameter \p in (22) depends on the direction 
of the incident wave and on the direction of the point 
of observation such, that \p (and consequently t ) 
tends to zero when one of the boundaries of geo- 
metric optics is approached. The asymptotic ex- 
pansion of the integral (21) for large yp can then be 



obtained by means of (18), (19), and (20): 

/(^)«Sc n r(n+i)( T p)- n -* 

71 = 

-fc- T#eo '<**> Erfc[(2 7 p)* sin *\ . (25) 

The coefficients c n are determined by (10) with 

g*(t)=g(t) + T' 1 2-h-^sin(^)(t+2sm^/2)-- 1 (26) 

and g(t) is given by (22). The function g*(t) fulfills 
the conditions (a) and (b) in the region — 7r<arg2<V. 
The asymptotic expansion (25) holds therefore in 
— 37r/2<arg (7p)<37r/2 and is valid also for ^=0 
(in which case all the c n vanish) and assumes in this 
case the correct value of (21) 



l(a,0)=-%e-^. 



(27) 



The results produced here are equivalent to those 
in [6,7]. 
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